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With time-dependent Lindblad operators, an open system may have a time-dependent 
decoherence-free subspace (t-DFS). In this paper, we define the t-DFS and present a necessary 
and sufficient condition for the t-DFS. Two examples are presented to illustrate the t-DFS, which 
show that this t-DFS is not trivial, when the dimension of the t-DFS varies. 
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I. INTRODUCTION 

In recent years, many proposals are presented to pro- 
tect quantum systems against decoherence. Except the 
method to weaken the coupling between the system and 
its surroundings, these proposals include the dynamical 
decoupling [IHH , quantum error-correcting codes [H, Q , 
the scheme based on the decoherence-free subspaces 
(DFSs) or noiseless subsystems (fUll, and the scheme 
based on the quantum reservoir engineering [ToL [Tl| . 

The DFS is defined as a subspace within which the 
system undergoes an unitary evolution [91 . I t was exper- 
imentally realized in a variety of systems [12h14| and has 
drawn much attention because its potential applications 
in quantum information processing [l^, [l6j]. Generally 
speaking, there are two ways to define the DFS. The first 
is given in Q, where the DFS includes all states that 
each state p satisfies C(p) = 0, where £(...) represents 
the Lindblad superoperator (we call it as the first defi- 
nition of DFS). The second definition was given in Q, 
which is formulated as follows. Let the time evolution 
of an open system with Hilbcrt space Hs be governed 
by the Markovian master equation. A decoherence-free 
subspace "Hdfs is defined as a subspace of Hs such that 
all states p{t) in DFS fulfill d t TT[p 2 (t)} = 0, for Wt > 
0, withTr[p 2 (0)] = 1. By this definition, it was proved 
that the subspace Wdfs = Span{|$i), |$ 2 ), • ■ ■, |$m)} 
is a DFS if and only if each basis of "Hdfs satisfies 



,M;a = 

2 ^-<a 



1, K, and "Hdfs is 
c* a F a -c a Fl) , where 



invariant under H c r = H 

H is the Hamiltonian of the open quantum system 

Notice that the basis of the aforementioned DFS is 
time- independent, the DFS is then time-independent. 
For time-dependent Lindblad operators, however, a time- 
independent DFS may not exist, then a natural ques- 
tion arises:What is the DFS for open systems with time- 
dependent Lindblad operators? If the DFS is time- 
dependent, what is the condition for the system to re- 
main in this subspace? 

In this paper, we will give a detail analysis for a time- 
dependent DFS (t-DFS). The analysis is given based on 



the Lindblad master equation with time-dependent Lind- 
blad operators. We shall adopt the second definition for 
the t-DFS and develop a theorem to give a necessary and 
sufficient condition for the t-DFS. 

This paper is organized as follows. In Sec [III we shall 
give the definition for the t-DFS and derive a necessary 
and sufficient condition for it. In Sec JIIIl we present two 
examples to illustrate the t-DFS, showing that the neces- 
sary and sufficient condition can be satisfied by manipu- 
lating the Hamiltonian. Finally we conclude by summa- 
rizing our results in Sec lIVI 



II. NECESSARY AND SUFFICIENT 
CONDITION FOR TIME-DEPENDENT DFS 

In this section, we will present a necessary and suffi- 
cient condition for the t-DFS and show that this condi- 
tion can be satisfied by manipulating the Hamiltonian of 
the open system. 

Consider a system S with ./V-dimensional Hilbert space 
"Hs coupling with an environment. The time evolution 
of the open system is assumed to be governed by, 



P {t) = -i[H(t),p(t)]+C(t)p(t), 

1 



F a {t)p{t)Fl{t) - -{Fl(t)F a (t),p(t)} 



(1) 



where F a (t) (a = 1, 2, 3..., K) are time-dependent Lind- 
blad operators and £(t) describes the Lindblad superop- 
erator of the open quantum system. 

Theorem 1. Let the time evolution of an open quan- 
tum system in a finite- dimensional Hilbert space be gov- 
erned by Eq.Uty with time- dependent Hamiltonian H{t) 
and time- dependent Lindblad operators F a (t). The sub- 
space 

H DFS (t) = Span^xit)}, |* a (t)>, • ■ •, \* M (t))} (2) 

is a t-DFS if and only if each basis vector of H.DFs{t) 
satisfies 



F a (t)\* j (ty=c a (t)\* j (t)),j = l,...,M;a=l, 
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(3) 
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and 'HDFs(t) is invariant under 
H eff (t) = G(t)+H(t) 

+U2{c* a (t)F a (t)-c a (t)Fl(t)). (4) 

a 

Here G(t) = iW (t)V (t) and U(t) is an unitary operator 

M N-M 

U(t) = '£\$ j (0))(<Z> j (t)\ + £ l^(0))<^(*)|. (5) 

j=l n=l 

Proof. Firstly, notice that the effect of U (t) is to map 
a set of time-dependent bases of Hs into a time- 
independent one. Transforming the density matrix p(t) 
into a rotating frame, i.e., p(t) = U(t)p(t)U 1 ' (t), we write 
the master equation Eq.{T]) as, 



am 
at 



= -i[H(t),p(t)] 

-iG{t)p(t) + ip{t)G{t) + C{t)p{t), (6) 



where, 



H(t) = U(t)H(t)UHt), 
mm = J2(F a (t)p(t)Fl(t) - l/2{F^t)F a (t),p(t)}) 

a 

with F a (t) = U(t)F a {t)W(t). Clearly, G(t) = 
U(t)G(t)W(t) = iU(t)W(t) = -iU(t)U^(t) = Gt(t), 
this indicates that G(t) is a Hermitian operator. By 
defining new Lindlbad operator as F a (t) = F a (t) — c a (t), 
the decoherence terms in Eq.(JB]) can be rewritten as 

C(t)p(t) = C(t)p(t) 

^(c*(^ a (t)-c a (t)Ft(t)),p(t) 



(7) 



where 

mm = j2( p mp(t)Fi(t) i/2{Fi(t)F a (t),p(t)})- 

a 

Note that the requirement of invariance of the sub- 
space %dfs(£) under the operator H e ff(t) in Theorem 1 
implies (^(0)|^e ff (t)|^(0))- (^(t)\H eS (t)\^(t))_= 
for Vn, j. In the rotating frame, the subspace "Hdfs 
spanned by {|<Fj(0))} must be invariant under the oper- 
ator 



B eS (t) 



H(t) + G(t) 
i 

+ 2 



l -Y^(<(t)F<*(t)-c<*(t)Fl(t)). (8) 



Now we prove that the condition is sufficient. Any 
state \ f(t)) e rlr>Fs{t) can be expanded by {\&j(t))}, 



M 



b(f)> = £«i(*)l*i(*)>- 



(9) 



By defining \ip(t)) = XJ(t)\ip(t)), p(t) and will be 

written as, 

m = u(t)\<p®){<p(WHt), 

C(t)p(t) = -i [i E a K(t)F a (t) - c a (t)Fl(t)),p(t)] , 

(10) 

where F a (t)\p(t)) = 0x \p(t)) = and C{t)p(t) = have 
been used. Hence the evolution of p(t) is governed by 



dp(t) 
dt 



-i[%s(t),p(t)] 



(11) 



|Tr[p 2 (t)] = §-Tr[p 2 (t)} = 2Tr[p(i)A(*)] 

= -2iTr{[R eS (t),p(t)]p(t)}=0- (12) 



Thus the conditions Eqs.(|3]) and ((4]) arc sufficient for 
JLTr[p 2 (*)]=0. 

In order to prove that the conditions are necessary, we 
assume that the set of basis {!<&;(£))} spans a t-DFS. At 
a fixed instant to, the quantum state \ip(to)} embeds in 
t-DFS 



lv(*o)> = E^ 1 Oi(*b)l*i(«o)>. 
P (t Q ) = |v(*o))M*o)|, 



(13) 



which implies 

= ^Tr[p 2 (i)]| t=t0 = 2Tr[p(t )C(t )p(t )} 

= 2(<p(t )\jC(t )p(t )\<p(t )). (14) 

Then, by using Eq.([T]), we obtain 

= (v(to)\£(t )p(to)Mt a )) 

= ]T 7q [(^(to)|F a (t)|^(io))(^(to)|Ft( i )|^(t )) 

(X 

-{tp{t )\Fi{t)F a {t)\y{t ))] . (15) 

Without loss of generality, we set F a (to)\(fi(to)) = 
c a {to)\f(to)) + \<P X (to)) witn l^^o)) being a (non- 
normalized) state orthogonal to state |<p(fo))- Submitting 
this into Eq. (jT5j) . we have 



(16) 



Since j a > for Va, we have (tp- 1 (to)]^ (to)) = 0. It 
straightforwardly follows F a (to)\tp(t )) = c a (t )\ip(to)} ■ 

Next, we prove by contradiction that the basis vec- 
tors {\<&j(to))} are the eigenstates of the Lindblad oper- 
ators F a (to) with the same eigenvalues c a (to). Suppose 
that two arbitrary eigenvectors of the Lindblad operator 
|$fc(t)) and |$fc/(t)) have different eigenvalues, i.e., 



F a (t)\$ k (t))=c atk (t)\$ k {t)), 

F a (t)|* fc /(t)) = Ca,tf(*)l*k'(*)>- 



(17) 
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The state \<j)(t)) = (|$ fc (t)) + |$ fc - (t)))/ V% must be not an 
eigenstate of F a (t). However, since \<j)(t)) is a state lying 
in %dfs(*), the state should fulfill {C[\<j>(t)) {4>{t)\}) = 0. 
Hence the eigenvalues must be equal for all basis vectors. 

On the other hand, according to Eq.([7]), for p(t) = 
U(t)\<p(t)){(p(t)\Ul(t) with F a (t)\ip(t)) = c a {t)\(p{t)), we 
have 



where T is the time-order operator. If \ip(t)) = 
H e g(t)\(p(t)) is still in subspace HdfS: |^(*)) wm always 
bein'HDFS- Thus, \(p(t)) will evolve unitarily in "Hdfs(£)> 
if H c ff(t)\tp(t)) is still a superposition of the basis vectors 
{|$j(t))} of the t-DFS. Hence the conditions are neces- 
sary for ^Tr[p 2 (t)] = 0. □ 

Remark. Examining the unitary transformation Eq.([5|), 
one may suspect that the t-DFS is the same as (time- 
independent) DFS, in the sense that a state in the t-DFS 
T^dfs (t) can be obtained from the time- independent DFS 
%dfs by the unitary transformation U(t). This is not 
true, because for a t-DFS, its dimension can change with 
time. When the dimension of the t-DFS changes, we can 
not transfer the time-independent DFS into a t-DFS by 
an unitary transformation. In Sec lIIIl we will present an 
example to illustrate this remark. 

We now show how to realize a t-DFS by the Theorem 
1. Suppose that there is a set of degenerated eigenstates 
of all Lindblad operators F a (t). As stated 
above, if "HdfsW i s a t-DFS, it must be invariant un- 
der H cS (t), i.e., 

($^(t)|H eff (*)|$ i (*))=0 J Vn,i. (20) 

Substituting Eq.flU) into Eq. ([2TJ| and considering 
F a (t)\$j(t)) = c a (t)\^j(t)), we obtain 

-l^ a c a (t){$i(t)\Ft a (t)\$M = 0. (21) 

a 

Once the evolution of the Lindblad operators F a (t) are 
fixed, the basis vectors of Hr>Fs(t) are specified. The 
task to achieve an unitary evolution in the t-DFS relies 
entirely on the Hamiltonian of the open system, 

- l ^^<^(t)\F a (mi(t)). (22) 



+C(t)p(t) 

= -i[H eS (t),p(t)}. (18) 



(19) 



In the next section, we will present two examples to show 
that the quantum state of the open system can be re- 
stricted to evolve unitarily in the t-DFS by manipulating 
the Hamiltonian according to Eq. (|2"2")l . 

III. EXAMPLES 

In this section, we will present two examples. In the 
first example, we show that the sufficient and neces- 
sary condition for the t-DFS provides us a way to keep 
the open system in the t-DFS by manipulating the sys - 
tem Hamiltonian. In contrast to the scheme in Ref . [17l] . 
where the adiabatic condition is required to maintain the 
system with high fidelity in the t-DFS, here the adia- 
batic condition is removed. Here we should mention that 
some effort has been done on such issue, in which the 
time-dependent Lindblad operators evolves unitarily [l8j]. 
However, for the t-DFS, such a limit on the evolution 
of Lindblad operators is not necessary, therefore t-DFS 
we proposed is general. From the other aspect, by the 
first example, we explain why the adiabatic condition is 
required in Ref.[l7j-the Hamiltonian in that paper does 
not satisfy the condition of t-DFS. It seems that the t- 
DFS can be found by an unitary transformation, this is 
not the case when the dimension of the t-DFS changes, 
this point will be shown in the second example. 

A. Driven E-type three-level atom coupled to a 
broadband time-dependent squeezing vacuum 
reservoir 

Consider a S-type atom coupled to a time-dependent 
broadband squeezed vacuum field 0, [2(| ■ For squeezed 
vacuum field, the initial state can be written as 

\vac(r/)) = K(i])\vac), 

where K(tt) is a multi-mode squeezing transformation. 
There are three classical fields ilj(t)(j — 1, 2, 3) interact- 



dt 



-iG(t)p(t) + ip(t)G(t)-i[R(t),p(t)] 



and 



p(t) = Tcxp 



i / H off (r)dT 



p(0)Texp 



i / -ff c ff(r)dr 



4 



n, 



i> 



|0) 



l-i) 



FIG. 1: (Color online) Schematic energy levels. A H-typc 
three- level system driven by three classical fields Slj(t). We 
treat this system as an open system since we will consider it 
coupled to a broadband squeezed vacuum reservoir. 



ing with the atom as shown in FIGQ] Under the Born- 
Markovian approximation, the evolution of the atom is 
governed by the following master equation , 
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-i[H(t),p(t)}-j{Rl(r))R(r))p(t) 
-p{t)Rl(rj)R(Tj) - 2R( V )p(t)R^r,)}, (23) 



where H(t) = [fi 1 (t)|l)(0|+O 2 (t)|0)(-l|+O 3 (t)|l)(-l| + 
h.c], and 



R(rj) = S'cosh(r) + exp(i<j))S^ sinh(r) 



(24) 



l*2"(*)) = I )- According to Ea.([2"2"]). the classical fields 
f2j(i) must be modulated as follow, 

Qi(t) = cosh(r) exp(iwoi), 
fl 2 (t) = sinh(r), 

S7 3 (t) = wo sinh(r) cosh(r) exp(iw i), (26) 

such that |$df(<)) forms a t-DFS. 

The numerical results are presented in FIGJ2J in which 
the initial state is chosen to be I'&dfIO)) = c(r)| — 1) — 
s(r)\l) with r = 1. We plot the purity P{t) = Tr[p 2 (i)] 
with the parameter uio = O.I7 (FIGUJa)) and uio = IO7 
(FIGHfb)) as a function of time. From FIG|3] we can 
see that the control fields fij(t) play an important role 
in this scheme. Without the control fields (solid line in 
FIGH]), the decoherence-free evolution can be approxi- 
mately realized only under the adiabatic limit, i.e., the 
smaller the ujq is, the better the state remains in the 
t-DFS. This can be understood by examining Eq.(f2"2"j): 
When the control fields are absent, the t-DFS condition 
is equivalent to that the right side of Eq.([22"j) is zero, i.e., 
ujq — > 0. In other words, to make sure that the evolu- 
tion of quantum state is decoherence free, the squeezing 
phase has to change adiabatically. As expected, when the 
open system are controlled by the classical fields flj (t) as 
Eg. (f2l)]) (dash line in FIG 12]), the evolution of the quan- 
tum state is always unitary regardless of how fast the 
squeezing phase changes. 



The operator S = |1)(0| + |0)(— 1| denotes the absorp- 
tion of an excitation from time-dependent broadband 
squeezed vacuum field, and rj(t) = rexp(icj>) is time- 
dependent squeezing parameter with polar coordinates 
<t> £ [0, 27r] and r > called phase and amplitude of the 
squeezing, respectively. Form Eq. ([2~4]) . the state 



|*dfM)) = c(r)|-l)-exp(i0)s(r)|l) 



(25) 



with c(r) = cosh(r) / \J cosh(2r) and s(r) = 

sinh(r) / y / cosh(2r) , satisfies i?(ry)|$DF(i)) = 0. Here 
we assume that the phase of squeezing parameter is 
time-dependent with a linear trend <p = u>ot, and the 
amplitude of squeezing parameter is constant. 

It has been shown in Refs.[ljl[2(| that the decoherence- 
free evolution can be achieved by adiabatically changing 
the squeezing phase. Comparing with Ref.[19(, we will 
show here that the adiabatic limit can be removed when 
the time-dependent Hamiltonian in Eq. (|23|) is manipu- 
lated according to Eq. (f22|) . In other words, the open sys- 
tem could evolve unitarily in the t-DFS by employing the 
scheme proposed in Sec [II] The master equation Eq. ((2"3"]) 
can also be realized in the model of a pair of trapped four- 
level atoms [2lJ, hence there are many manner to realize 
it in experiment. 

In the following, the state |$df(£)) is chosen as the 
basis of the t-DFS, whose complemental space is then 
spanned by |$i"(f)) = s(r)\ — 1) + exp(iu>oi)c(r)\i) and 



B. A toy model for DFS with time-dependent 
dimension 

In this subsection, we present a toy model with vary- 
ing dimension of t-DFS. Consider a five-level system cou- 
pled to two different broadband squeezing vacuum fields 
and six classical control fields (f2j and f^), as shown 
in FIG[3] a(ui) and a'{ui) are the annihilation operators 
of the environment with different polarization and fre- 
quency lu. Suppose the squeezing vacuum states for the 
modes a(tu) and a'(tu) have different squeezing param- 
eters r/i(t) = r\ exp(jwoi) and 772 (t) = r 2 exp(iwot), re- 
spectively. Under the same assumptions used in Eq. 
the master equation can be written as (l9j . 



P(t) 



-i[H(t),p(t)} 



F a (t)p(t)Fl(t) - -{Fl(t)F a (t),p(t)} 



(27) 



where the Lindblad operators are F a (t) = cosh(r Q )S' Q + 
exp(wot) sinh(r ce )S't (a = 1,2) with Si = |1)(0| + |0)(-1| 
and 5 2 = |1')(0| + |0)(-1'|. The Hilbert space Hs(t) of 
open system can be spanned by the following orthogonal 
normalized bases, 



|$dfi(*)) = ci\ - 1) - exp(iw t)si|l), 



(28) 
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FIG. 2: (Color online) The purity P(t) as a function of coot (in 
units of 7r) with time-dependent Hamiltonian (dash line) and 
time-independent Hamiltonian (solid line) for (a) ujq — O.I7 
and (b) uio = IO7. The the other parameters chosen are 7 = 1, 
r = 1. 




ii> fl'i/ 

/ / a'(co) 



ft 2 



a'(a)) 



i- i-i') 



FIG. 3: (Color online) Schematic energy diagram. A five-level 
system, transition |1) -H- |0) and |0) <-> \ — 1) are driven by 
mode a, |1') <-> |0) and |0) o | — l') are driven by mode a', 
and |1') <-> I - l'> by 6. 



|$DF2(*)) = c 2 
l*ft*)> = *l| " 
1*2 (*)> = *a| - 
l*o"(*)> = |0>, 



-1') -exp(zwot)s 2 |l') 
1) + exp(iw t)ci|l), 
l') + exp(iLu t)c 2 \l') , 



(29) 
(30) 
(31) 
(32) 



with Si=sinh(ri)/ yj 'cosh(2rj) and 

Cj= cosh(ri)/ y/ cosh(2rj) . It is not difficult to check that 
this model admits a two-dimensional t-DFS at most, 



FlWI^DFl) 
i ? l(t)|*DP2) 

The Hamiltonian 



F 2 (t)|$DFl> 
-F2(*)|*DF2> 



0. 
0. 



(33) 
(34) 



H(t) = H (t) + T(t)fi(|*DFi(t))($DFa(t)| + /i-c), 

(35) 

includes two terms: The first term H$(t) is required to 
construct the t-DFS with time-dependent dimension, 

Ho(t) = (Oi(i)ll><0| + O 2 (i)10)<-1[ + 3 (t)|l><-l| 

+n 4 (t)|i')(o| + n 6 (t)|o)(-i'| + n 6 (t)|i')(-i'| 

+h.c), (36) 

where the classical field strengthes, according to Eq.(f2"2"j). 
are designed as follow, 

f2i(i) = cosh(rx) exp(iu>ot), 

02(i) = sinh(ri), 

f2 3 (t) = uj sinh(ri) cosh(rx) exp(— iuiot), 

Fl'i(t) = cosh(r 2 ) cxp(iui t), 

( 0, LJ t <tt/2; 

n' 2 (t) = </ (1 + cos(2u;oi))sinh(r 2 ), tt/2 < oj t < n; 
[ sinh(r 2 ), uj a t > n. 

£l' 3 (t) = too sinh(r 2 ) cosh(r 2 ) exp(— iujQt). (37) 

This implies that, for 0<uj t<Tr, the t-DFS %dfs(*) 
is one-dimensional and spanned by |$dfi(^)), since 

(^(t)|H eff (*)|$ DP1 (*)) = ($DF2(t)|H e ff(t)|$DFlW> = 

0. For w Q t > 7T, the DFS %dfs(*) is two- 
dimensional spanned by |$DFi(i)), and \^BF2{t)}, since 
($t(t)\H eff (t)\$ DF1 {t)) = ($ ? ± (t)|H cff (t)|$ DF2 (t)) = 0. 
Thus, the dimension of the t-DFS "Hdfs(^) changes with 
time. To illustrate the dimension changing with time, 
we introduce the second term of the Hamiltonian that 
induces a transition between |$dfi(*)) and |$df2(£))> 
where T(t) = 9(ujQt — n) is a step function and fl is 
transition coefficient. This means that the transition is 
allowed for Lo$t > 7r but it is forbidden for < toot < n. 

The numerical simulations are presented in FIGJ3J The 
population on |$dfi(0) (dash line) and |$DF2(£))( son d 
line) , 



Pl(t) = (*DFl(t)|p(t)|*DFl(«)>, 
P2(t) = <*DF2(*)|/0(*)|*DF2W>, 



(38) 



are shown in FIG|2|a). As predicted, when < ujot < ir, 
the population stays on |$dfi(£)); when u^t > 71", the 
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FIG. 4: (Color online) (a) The evolution of population on 
|$DFi(t)} (solid line), |<3?DF2(i)} (dash line) and population in 
%dfs(£) (dots line) with n = r% = 1 and f) = wo = 7 = 1. 
(b) The evolution of the purity P(t) with T(t) = 1 (dash line) 
and Tit) = 6(u) t - n) (solid line). 



population transits between |3>dfi(£)) and | < £ , DF2(i)) i but 
the total population Pi(t) + P 2 (t) (dot line in FIGHfa)) 
remains unchanged, this means that the system does not 
leak out the t-DFS. If the Hamiltonian does not satisfy 
the t-DFS condition, for example, T(t) = 1 instead of 
T(i) = 9(u)ot — ir), the purity P(t) will decay (dash line in 
FIGJlfb)). This implies that the dimension of the t-DFS 
really changes in the time evolution at time to given by 
ujoto = 7r. Once the dimension of the t-DFS changes, the 



time-independent DFS and t-DFS can not be connected 
by an unitary transformation. In the other words, the 
t-DFS is not trivial. 



IV. CONCLUSION 

In summary, we have defined and presented a nec- 
essary and sufficient condition for the time-dependent 
decoherence-free subspace (t-DFS) for open systems. In 
contrast to the time-independent DFS, the basis of this 
t-DFS is time-dependent. Besides, the dimension of t- 
DFS may change with time, this implies that we can 
not trivially get the t-DFS by unitary transformations. 
Two examples arc presented to illustrate the t-DFS. In 
the first example, we show in details how to manipu- 
late the Hamiltonian of a S— type system to realize a 
one-dimensional t-DFS, while in the second example, we 
through a toy model to show that the dimension of the 
t-DFS can change with time. The later indicates that the 
t-DFS can not be derived by an unitary transformation 
from the conventional DFS. 

The observation of the t-DFS and the prediction made 
for the t-DFS is within reach of recent technology. In fact, 
in a recent proposal [22j , the authors proposed a scheme 
to entangle two atomic ensemble of Cesium at room tem- 
perature by engineering reservoir |23 |. 24j. These tech- 
niques together with measurement |23(, can realize the 
t-DFS. For example, the time-dependent Lindblad oper- 
ators may be achieved by modulating the detuning be- 
tween the pumping field and the atom or by modifying 
the Zeeman splitting. (A detail relation between Lind- 
blad operators and parameters of environment can be 
found in Refs. [22j.) The parameters in the Lindblad op- 
erators can be modulated in experiment (25|. All these 
together can lead to the t-DFS. In practice, the Lind- 
blad operators are not known a priori. However, we can 
re-design a decoherence-free subspace by controlling the 
Hamiltonian and engineering the reservoir. This means 
that we can improve the decoherence by designing a DFS 
to include the states of interests. 
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